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NOTATION 


The following syitibols are used in this paper: 
b Bottom 

G Concentration of Species in the Water Column, M/L^ 

D Depth of Water in the Bay, L ‘ 
fi Rate of Mass Transfer by Evaporation, L/T 
Dispersion Coefficient, L^/T 
f Bay Bottom Friction Factor 
g Gravitational Acceleration, L/T^ 

H Height of Water above a Cell Datum, L 
K Constants 
o Source or Sink Term 

Q Discharge Rate, L^/T 

R Rate of Mass Transfer by Rainfall, L/T 

Rate of Disappearance or Appearance of Mass, M/L^T 
r Resuspension 

s Settling 

Surface 
t Time, T 

V Resultant of the Velocity Vector, L/T 

V Local Grid Velocity, L/T 

W ' Angular Velocity of Earth, L/T 
X Distance (East-West) , L 

y Distance (North-South) , L 

z Wind Speed, L/T 
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NOTATION (Continued) 


n Wind Speed, L/T 

0 Temperature 

cp Wind Direction 

Angle Measurement in the Coriolis Term 

Differential Operation 
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ABSTRACT 


There is a growing awareness that the natural resources of the 
world are limited. This fact gives credence and a sense of immediacy to 
those who are trying to better understand and describe those processes 
which affect the amount and the quality of these resources. One of the 
most abundant and most taken-for-granted natural gifts is water. The 
waters adjacent to coastal regions are some of the most often studied 
because of their importance to man. 

The coastal environment is a vital part of man's daily activity - 
providing food, recreation, jobs and habitats. Thus the already com- 
plex, dynamic natural processes which maintain a balance between fresh 
water and saline water is further confounded by man-made impacts. To 
minimize adverse events on these areas, a clear understanding of the 
properties and behavior of these systems must be established. Plans 
formulated with technically sound data are far more likely to produce 
results which are both environmentally and economically sound. 

In recent years, studies have been accelerated to better charac- 
terize the coastal waters and to better describe the processes which take 
place in these areas. Studies have included models - both mathematical 
and physical - as well as the more traditional investigations involving 
data acquisition - both field oriented and remotely sensed. The 
interaction of these methods provide techniques for the rapid prediction 
of changes in the system and the impact that these changes have on 
water quality and behavior. 

The National Aeronautics and Space Administration, Marshall Space 
Flight Center (MSFC) has had a continuing interest, through their 
Environmental Applications Branch, in the application of remote sensing, 
automatic data processing, modeling and other aerospace related 
technologies to hydrological engineering and water resource management. 
One such investigation includes; major hydrological and water resource 
emphasis for the entire river drainage system which feeds the Mobile 
Bay estuary. This study has a particularly significant and timely 
purpose as a result of the development of the Tennessee-Tombigbee 
Canal, a project which will connect the Tennessee-Ohio River systems 
to Mobile Bay via the Alabama-Tombigbee Rivers. The impacts created by 
this construction and other such developments on the water quality and 
hydrological characteristics of the bay are of primary concern and 
importance . : - - 

As part of that study, MSFC has funded projects under contract 
NAS8-29100 to investigate the adaptation and implementation of existing 
mathematical modeling methods for the purpose of describing the be- 
havior of Mobile Bay, Of particular importance are the interactions 
that system variables such as river flow rate, wind direction and speed, 
and tidal state have on the water movement and quality within the bay 
system. 
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Specific achievements of these modeling studies include; 

1. The development of a rapid, predictive technique for 
establishing baseline conditions within the bay system 
from which other studies can be compared. 

2. The determination of the dynamic interchange occurring 
between the bay and rivers, the bay and the Gulf of 
Mexico and Mississippi Sound, and the bay with 
adjacent land masses. 

3. The investigation of bay physio-chemical processes 

and the impact that these have on the water mass. Such 
items as changes in river flow rates, runoff and 
pollution loading are assessed in terms of the way 
they affect distribution and transportation properties 
■within Mobile Bay (currents and tidal elevations) . 

4. The determination of material transport behavior within 
the water system as controlled by system variables 
(salinity, sediment transport, coliform bacteria) . 

5. The establishment of a basic model useful to extend 
present capabilities to include other material trans- 
port problems (BOD/DO distributions, oil spill transport). 

6. The development of a method to interface the bay model 
with selected subsystems designed to provide a closeup 
view of certain local disturbances (dredge discharge 
material transport. Little Dauphin Bay model). 

7. The establishment of a method with the capability of 
interacting with field oriented, data processing and 
remote sensing programs which are concurrent with 
contract NAS8-29100. 

These achievements and the programs which are suggested for further 
investigation are the subject of this final report. The work presented 
represents nearly three and one half years of investigation at The 
University of Alabama under the support of contract NAS8-29100 with the 
Marshall Space Flight Center of the National Aeronautics and Space 
Administration. 
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The primary objactive of ocntract NAS8-29100 was the demonstration 
of thi ntlllty of mathematical modeling methods i 

1, the rapid, predictive 

in a complex system such as Mobxle Bay, 

1 nf a strong data base from which the 

2. the development of a § auality and movement 

impact of future events on bay quatxry a 

could be compared, and, 

3 the establishment of a technique to 

sensS and automatic data processing methods. 

TO achieve these ohiectives “f ^,,7” :rifn«1fd 

out the course of the study, leted elements. Each of these 

Stiheif rtatur:n?thrircontrlbution to achieve the contract 
objective are listed in Table 1. 

nf t-he report, detailed results derived from 
In subsequent sectxons of the repo , presented to illustrate 

the modeling studies (listed xn PriL to that, however, a 

the range, depth and utx_x y system, a summary of the model 

brief description of the ® jata used to Implement the 

aquations, and the sources Investigations into a 

models are needed to place cne _„view. 

familiar setting for dxscussxon and revxew. 
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Table 1. -Contract Plan by Work Element Shoving Current Statue and Past Reports leeued. 


Work Element 
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1974 

1975 
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B 

B 
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B 

B 

B 

B 

B 
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B 

B 

!B 

!BI 

I. Review of Existing Mathematical 

Models for Estuarine Systems 

1 

1 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■1 

11. Adstptsrtlon of Models to Ifaiblle 
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1X1. Data Acquisition/Identification 
for Calibration, Verification, 
Parameter Seoeitiylty Analyses 

II 
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1 
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1 

1 

1 

1 
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1 
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IV. 


L . 


V. 


c. 


ente 


Completed under ikrlne Science 
Programs Project. 


II. 


II. B. 3, 


One third completed under Marine 
Science Programs Project; two 
thirds under contract MRS8-29100. 
Model developed for BOD-DO; lnsuf< 
flclent data for verlflcatlon/cal' 
ibratlon. 


D. Conform Bacteria 


g, 


IIIJl. 


III.B. 

J III.C. 
IIl.D. 

III.E. 


Application of Developed Models 
to Specific Water Quality and/or 
Bay System ProblecM 

A. River Discharge, Wind, Tidal 
Conditions (Wormal) 

B. Pollution Loading fCollform) 


C. River Discharge, Wind, Tidal 

Conditions (Abnormal) 

D. Subayatem Modeling 
1. Little Daaphln Bay 

Bay Dradi^ Dis cha r g e 
3. Material TraBeportCeerir 


A. Baaotti 


Tr a w efe r 




Continuing effort due to Interact- 
ion of hydrodynamic model and 
isaterlal transport sndels. 

Tidal cycle average data only 
(literature) , 

Insufficient data. 

Monthly average data only (State 
Depertment of Health). 

Seasonal and longer periods (llt- 
eratnre) 


Abnorsul (extreme river flows, 
storm surges and winds); Insuffi- 
cient data for verlflcatlon/call- 
bratlon. 

Little Dauphin Bay; Insufficient 
data for verlftcatlon/callbration, 
Maintenance Dredging activities; 
Corps of Engineers data. 
Incorporated to Increase users 
participation (replaced elements 
III.C. and IV.C.). 


Reports Issued 


BER 168-112 

BER 169-112 
BER 174-112 
BER 185-112 
Rgg 2 na.n? 


KR 168-112 
BER 169-112 

BER 185-112 
BER 208-112 


interim reports (fourth 
1976 Includes a final report). 


BER 168-112 
BER 169-112 
BER 185-112 


BER 203-112 

BER 208-112 
BER 203-112 


BER 209-112 
(Pinal) 
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THE MOBILE BAY SYSTEM 


Mobile Bay is approximately 400 square miles in area and is located 
on the northeastern shoreline of the Gulf of Mexico east of the 
Mississippi River delta. The estuary is about 31 miles long and varies 
in width from 8 to 10 miles in the northern half to about 24 miles wide 
in the southern portion (Figure 1) . The southeastern region of the 
estuary is referred to as Bon Secour Bay. The southern end is blocked 
from the open Gulf by land barriers; Gulf Shores to the east and Dauphin 
Island to the west. There are two passes located in this area also; 
the main pass which connects with the Gulf at Mobile Point and the pass 
which connects with Mississippi Sound at Cedar Point. 

The bay is the terminus of the Mobile River system which consists 
of more than 43,000 square miles of drainage basin; the fourth largest 
in the United States. Variations in river discharge rate from a ten 
year average high of 136,000 cubic feet per second in March to a ten 
year average low of 16,000 cubic feet per second in September have been 
recorded. During periods of heavy runoff the bay receiving waters are 
dominated by the high river flows to the extent that salinity intrusion 
within the bay is suppressed to the mid and lower reaches of the estuary 
Conversely, during low flow periods, salinity concentrations of 15 to 
20 ppt are measured in the upper bay. 

Also, as a result of the large region drained and the relatively 
high river flows during the rainy season, suspended sediment loads 
equivalent to about 5.5 x 10® tons/yr are carried to the bay. Estimates 
show tuat approximately two-thirds of the load is deposited within the 
bay resulting in long term bathymetric modifications controlled by the 
hydrodynamic and material transport properties of the estuary. The 
variation in sediment load observed during the period 1952-1963 is from 
a high of 1.0 x 3.0® tons for the month of March to a low of 5.5 x 10^ 
tons for the month of August (1) . 

Shoaling in the bay has averaged about 2 feet per century. 

However, there are portions of the bay which are highly stable and other 
regions which have rates of nearly 10 feet per century. These wide 
variations are a result of the complex, natural circulation patterns and 
man-made influences such as channel construction and maintenance dredg- 
ing activities which exist in the bay. Nearly 1.8 x 10® tons of 
suspended sediment bypass the bay annually. This material discharges 
into the Gulf through the two natural passes in the southwestern area 
of the estuary. During tidal dominant periods (especially during flood 
tide cycles) solids can be introduced into the bay from the Gulf. These 
materials are transported to the bay mouth by the predominantly east to 
west littoral current which occurs in the northeastern Gulf of Mexico. 

In addition to the riverine and tidal influences on circulation and 
material transport in the bay, wind direction and speed is also an 
important variable (2) . 
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Sewage, industrial waste disposals, and storm water overflow dis- 
charged into Mobile River and surrounding creeks from the Mobile 
metropolitan area, and excessive concentrations of bacteria in the 
Mobile River, result in the pollution of Mobile Bay. One method for 
expressing the bacterial content of these waters is to determine the 
total conform bacteria group count which gives an indication of the 
disease carrying bacteria or pathogenic content in the water. Because 
of this pollution, Alabama, under state laws and the regulations of the 
State Board of Health, periodically closes the bay to oyster harvesting 
as a safeguard to human health. 

Rapid predictive methods supplemented with spot analytical support 
could result in substantial savings of time and effort in analyzing bay 
behavior in the above categories. The method could also provide answers 
related to the abatement and prevention of serious upsets to the system. 
This study provides such a method which has a basis the application of 
conservation of mass and species equations subject to the bay ecosystem 
constraints. For this purpose, a two dimensional (surface), non- 
conservative species transport model is developed for Mobile Bay. The 
model is solved with a finite difference method and implemented by 
computer solution using a IINIVAC 1110 system. The hydrodynamic model 
for Mobile Bay developed by Hill and April (3) is used to provide basic 
current and dispersion coefficient data required by the non-conservative 
species transport model (NCSTM) and the conservative species transport 
model (CSTM) . Calibration and verification of the CSTM for salinity 
and sediment distributions, and the NCSTM for coliform bacteria were 
achieved using available field data for each species investigated. 
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MATHEMATICAL MODELS OF THE MOBILE BAY SYSTEM 


The equations which make up the mathematical modeling package for 
Mobile Bay are shown in Table 2. Modifications at the boundaries of 
the system, including river-marsh interaction, and stratification (salt 
wedge) effects, have been incorporated into the solution method. 
Similarly, coliform die-off rate constants and sediment resuspension 
and settling rates are introduced when these elements are studied with 
the model. The model results are presented as average concentration dis- 
tributions within the Bay corresponding to the time frame over which 
field data were available for calibration and verification. 


Verification Data Bases For The Mobile Bay Models 

Hydrodynamic Model - Synoptic hydrodynamic data at locations within the 
Mobile Bay system were received from the Corps of Engineers, Mobile, 
Alabama, for May 15 and 16, 1972. This information consists of tide 
charts and discharge rates experimentally determined over a thirty-four 
hour period. Tide heights are taken from the appropriate charts and 
converted to read from the model reference plane (mean sea level) . 

Fourier series are fit to data by the method of least squares. 

Manning coefficients vary from 0.010 to 0.050. A coefficient of 
0.050 is used in the marsh area to simulate the low flow rates expected 
in this zone. Values within the bay proper range from 0.010 to 0.018. 
Initially, efforts were made to account for variations in roughness 
created by oyster beds, channels and spoil banks. However, large changes 
in the Manning coefficient caused only minor changes in the flow on the 
scale of the model used. 

The hydrodynamic model is exercised over four tide cycles beginning 
with estimates from a previous run. The first step in the verifi- 
cation process is a check of the tidal heights at Mobile State Docks , 
Great Point Clear, Fowl River, and Bon Secour (Figure 2). Both tidal 
amplitudes and phases checked closely with the actual data. This is 
particularly significant in view of the fact that the forcing function 
in the Gulf of Mexico and at Cedar Point are smoothed data derived from 
storage equations as previously discussed in relation to boundary con- 
ditions. Other factors which may influence somewhat the exactness of 
the fit are localized winds and adjacent marsh areas that may flood at 
the high tide. 

The second verification step consists of a comparison of dis- 
charges at Main Pass and Cedar Point with field measurements taken by 
the Corps of Engineers (Figure 3) . Discharges are calculated by the 
Corps from periodic measurements at various locations in these passes 
at a depth of 0.2 and 0.8 times the depth of flow. An arithmetic 
average of these values is considered to be the average value for that 
location in the vertical direction. Horizontally, the area covered 
is half the distance to the adjacent measurement location on either side. 
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Table 2.— Mathematical Representation and Operational Modes of the 
Physical Models for Mobile Bay (3,9) 


Name 

Eouation Form 

Results 

Modes 

Continuity 

3Qx + ^ (R + E) 

3x oy St 

Tidal Height 

Tidal Cycle 
■Daily Avg 
Monthly Avg 
Seasonal 

Momentum 

x-Component 

+ gD ^ = KV^cosij' - fQQxD"2 

3t ox 

+ Qx(2Wsincp) 

X- Component of 
Surface Current 

Tidal Cycle 
Daily Avg 
Monthly Avg 
Seasonal 

y-Component 

1 

1 

+ gjj ^ fQQvD"^ 

3t oy ^ 

+ Qy(2Wsincp) 

! y-Component of 

Surface Current 

Tidal Cycle 
Daily Avg 
Monthly Avg 
Seasonal 

Species 

Continuity 

|E+%|e+Vy|£ -E(g + |^, 

Concentration of 
Species 



” ^(CVz(2s) - CvzCxj,)) 




+ R„ 



Salinity 

Ro = 0 

Salinity Concen- 
tration 

Daily Avg 
Seasonal 

Colrtorm 

Ro = KC; where K = f (0) 

Coliform Bacteria 
Concentration 

Monthly Avg 
Seasonal 

Sediment 

Ro - Klf(vs) + K2f(E) - K3f(v^) 

Suspended Sediment 
Concentration 

Seasonal 
Tidal Cycle 


Note: f in the above equations denotes a functional representation of the 
variable listed in parentheses. 





Slni f I'ocl'!: 


3.2 ~ 


b) Point; Clcni: 




o.£i~ 


• / v 

I / \\ 

'/ V 

'/ \ 

// \\ 


■ I ■ -o.fl-uj 

1230 1230 


_ c) Fowl River 


d) Bon Secour River 


0/|3() 203f) 

M.ilil'viv TilT'«.> 


_lJ_j -0. ll— J — •- 

1230 1230 


0630 2030 

Military Time 


Figure 2. 


Tidal Cycle at Four Loraciona in Mobile Bay for Model (Dashed Line) 
and Field (Solid I.ine) Results, 
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The correlation between actual data and model predicted data at 
Main Pass is excellent. This x<ras expected since flows through this pass 
are well behaved and influenced little by the flow at Cedar Point. There 
is a slight deviation at Cedar Point between model predicted values and 
field data. During the course of this study, it became evident that 
flows through Main Pass have a direct influence on the flows through 
Cedar Point (4) . 

A final step in the verification process is qualitative in nature 
and involves a visual comparison of flow patterns predicted by the model 
with those postulated from field observations (5 )• Comparative inflow 
(flood) and outflow (ebb) profiles are shorn in Figure 4. The -trend of 
water movement predicted by the model is in agreement with the literature 
data. Water entering through the Main Pass on the flood tide sweeps 
through Bon Secour Bay then turns to a northerly direction and flows 
toward the Blakely and Apalachee Rivers. During ebb flow, the water 
movement is from the upper bay and Bon Secour Bay toward the Main Pass. 
The patterns illustrated using literature data are derived from data 
collected over different periods of time. Therefore, they are only 
used here to indicate a trend in surface water movement for comparison 
xjith the depth averaged model results. 

Salinity Model - Considerable salt concentration data were found in the 
literature for the month of October, 1952. Average river flows obtained 
from the United States Geological Survey for this time period indicated 
a fresh water input of 12,000 cubic feet per second. The hydrodynamic 
model was exercised using the modified river flows to compute the 
pertinent data for the salinity model. 

Dispersion coefficients and net velocities, calculated in the hydro- 
dynamic model for each grid location, were used as input data for the 
salinity model. The Gulf boundary saline concentration was set at 
thirty-five parts per thousand. The concentration at Cedar Point for 
October, 1952, was elucidated from the literature (3) and set at twenty- 
five parts per thousand. Dog River and Mobile River were set at zero 
concentrations . 

The salinity wedge was accounted for in the Bon Secour area. This 
area was chosen rather than the ship channel for several reasons. 

First of all, data available indicates that the effects in the ship 
channel are minimal which may be due to the low comparative surface 
area involved on the scale of the model studied. Secondly, the 
literature (5) indicates that a large area in Bon Secour is influenced 
by the wedge. This is expected as a result of the flow patterns in the 
area. Finally, the model indicates that the salinity wedge in Bon 
Secour significantly contributes to the overall salinity patterns. This 
was achieved in the model using a first order equation for the rate of 
mass transfer from the salt wedge to upper water layers. Simulating 
the three dimensional salt wedge effect in this manner gave model results 
in closer agreement with reported field data (Figure 5) . 

The salinity mdoel was exercised for thirty-two tide cycles 
beginning with estimates from a previous run. Data from the 
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literature (5) were averaged for ebb and flood tides as well as in the 
vertical direction and compared with model trends. Results are in 
general agreement with the field data (Figure 5). A deviation along the 
stern side of the bay may be the result of unidentified fresh water 
flows in this area, as was Dog River at the outset of this study. This 
IS surmised from several pieces of information. Data from the 
iterature indicate a rather strong net outflow along the extreme 
western shore even for low fresh water flows (5) . iLs is in contrast 

Troflll Satellite photography and coliform 

P files which indicate the main thrust of net outflow is doxcm the ship 
channel and minimum flows adjacent to the western land boundary (6)! 

the possibility of additional fresh water flows, the model- 
predicted isohalines appear reasonable and lead to a study of natural 

SSiSoL?"''"'""' “ various^wind and river 


^ liform Bacteria Model - Total coliform group concentration data for 
various locations in Mobile Bay were collected by the Alabama State 
Department of Health for the period from January 1962 to August 1962 
Coliform group concentrations are obtained by analysis as described in 

Significance of EC Positive Organisms in 
Gulf Shellfish Growing Waters" (7) . 

The model is verified on a monthly basis, i.e. monthly average 
conditions are used, and the model results are tabulated and compared 

verlflc^rr criterion for model 

t-K XS based on how well model-predicted results fall within 

the field data range at the several locations within the bay for any 
given monthly period. ^ ^ 

f-K dependence of the species continuity equation on 

he hydrodynamic model of Mobile Bay for current distributions and 

coefficients, the first step in the verification procedure 
volves specification of data necessary for the proper description 
of the hydrodynamic behavior of the Bay. This includes the calculation 
mon y average river flow rates, wind conditions and tidal con- 
ditions for the period for which total coliform group concentration 
data are available (Table 3) . concentration 

thp total coliform dieoff rate constant used in 

the model is calculated as a function of monthly average water 

temperatures are estimated 
Bv R bimonthly average water temperatures of Mobile Bay compiled 
In Sn B ^ ^ recognized that water temperatures are not uniform 
in the bay. The degree of mixing that occurs between sea water and 
river _water within the bay will affect the temperature distribution. 

In this study temperatures are considered homogeneous throughout the 
ay. Temperatures can be adjusted linearly between the values corres- 
ponding to Gulf of Mexico water temperature and river water temperature 
to approximate real system behavior. In this study where monthly 
average values are Investigated the sea water intrusion effect can be 
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Total conform group concentration data for locations recognized as 
having severe pollutant input into the bay are used as loading con- 
c^^Katlons attach relevant grid cell (Table 5). They are held constant 
throughout each computation. Loading at the Mobile Rxver has been found 
to be the main source of pollution of Mobile Bay (9) . 

Results are presented as model calculated total coliform profiles 
within Mobile Bay (Figure 6 to Figure 13) . Similar results are tabulated 
for each month from January to August, 1962, during which the veri- 
fication phase is performed (Table 6 to Table 13). Total coliform con- 
centration vs. time (month) curves for the bay locations indicated in 
Figure 14 are also presented to show the trend in the concentration 
changes with season (Figure 15 to Figure 2?) • Details of the cali- 
bration and verification methods used can be found in reference (9) . 

Sediment Model - Mobile Bay experiences seasonal variaticDns in rainfall, 
runoff and sediment loading which can be broadly classified as low, 
medium and high in a manner outlined in Table 14. 

During the heavy runoff period (high) the bay receiving waters are 
dominated by the high river flows to the extent that salinity intrusion 
within the bay is suppressed to the mid and lower reaches of the 
estuary. This condition also causes the most severe material transport 
and deposition of sediment within the bay system. Conversely, during 
periods of reduced river flow and sediment loading, bay currents are 
dominated by the tidal influence. This results in a greater potential 
to transport sediment although the total volume is significantly reduced 
because of the low runoff conditions. 

The data used for this study are seasonal average 
sediment loads obtained from the literature for the period 1952-1963. 
This information is used in the hydrodynamic model to calculate a tidal 
cycle average velocity which is then related to sediment transpc^r 
potential by correlation with a critical velocity. Using this infor- 
mation, trend analyses of Mobile Bay sediment deposition reported as a 
potential function are included. These results include seasonal and 
longer variations in the sediment patterns . 


Table 5. Loadings of Total Coliform at Various Locations in MPN per 100 ml 






































•Figure 8. Model Calculated Total Coliform Concentration 
Profiles for March, 1962. 
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Table 7 . Total Coliform Concentration for Mobile Bay - 
February 1962 

Loading at Mobile River Mouth = 1/4 
TC 3 ^ (D.F. = 5 ) 

Correction Factor for E = 500 
= 0.29 day“^ 


Measured Data 


1 Station 
1 

No , 

i 

: 

i 

1 

1 

i 

No. cf 1 Monthly 

70^ Confidence Range 

Model 

Calculated 

1 

Sampling 

Average 

! X 

X - tSx 

X + tS^ 

Result 

i 

I -^-5 

2 

4,300 

338 

8,662 

1 , 981 


23,000 

18,125 

27,875 

15,415 ' 

1 

to 1 4 

5,000 

3,312 

1 

6 , 688 

1 

5,273 1 

6c ■ 4 

17,000 

10.750 

23,250 ! 9,677 I 

1 

: 6: i 4 

63,500 

41 ,000 

i 1 

86,000 i 10,548 I 

■■■ 

27.500 

^ r ^ ^ 

13,750 i 41,250 1 10,549 i 


881 i 2,419 1 

■3 L 

1 

66 : 4 1 . 8,000 

7,531 i 8,466 : 6,421 : 

. 

Cp \ h 

51.500 

33,562 I 69,436 8 , 40 ? ; 

' 75 , 1 ^ 

15,000 

1 

7,312 i 

22,686 1 5,302 ' 

I 83 1 3 

1.100 

1 

300 

1,900! 1 , 491 

1 1 i 

£S i ^ 1 

5,300 

! i ! 

1,800 1 8,800 . 3,783 i 

I 112 i 3 

1,380 

603 j 2,156 1 709 ' 
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Measured Data 


Station 
1 No. 

\ 

! 

■ 

f 

Nc. of 
Field 
Sampling 

Monthly 

Average 

X 

7O/0 Confidence Piange 

1 - 

Model ; s 

1 1 

Calculated 1 

fiesult j A 

■ k 


X + tSjj 

33 

2 

8,000 

4,126 

11,874 

3,938 ; 

36 

3 

25,000 

13.638 

36,362 

1 

80,734 ' 

59 

4 

1 

160,000 

-27,500 

347,500 

20,108 i 


4 

69,500 

49,688 

89,313 

41,863 I 

61 

4 

35,000 

15,000 

^ 

55,000 

43,519 ; 

62 j 4 

14,000 

i '■ "" 

6,625 1 21,3751 35,070 r 

63 i L 


r**” — i 

4,060 I 4,260 ! 11,815 P 

i 

66 1 4 

i 

36,000 

11 , 625 

60,375 1 

25,338 ;P 

67 , 6 

19,250 

■ 

8,625 

29,875 1 31,766 1 

75 i 

4 

15. 750 

1,763 

29,738'! 17,326 1 , 

1 

83 1 

4 

255 

186 

324 1 3,159 i 

I 

88 1 

4 

2,800 

1,363 

4,283 I 5,375 1 ^ 

« 1 

112 

3 i 

55 

-3 

113 1 

1,089 i 


































































Table 9. Total Coliform Concentration for Mobile Bay 
April 1962 

Loading at Mobile River Mouth = 1/4 
TC33^ (D.F. = 4) 

Correction Factor for E = 500 
da 


Station 

No. 


Measured Data 


No. of Monthly 70^ Confidence Range Calculated 
Field Average : Result 

Sampling x - tS^ 3c + tS^ 



1,340 


76,600 


162,000 


7,250 


27.500 


17,000 


7,100 


8,100 


2,750 


1,008 


2,072 


44,669 108,531 


-50,500 374,500 


5,688 


14,313 


7,000 


4,263 


2,725 


2,063 


8,813 


40 , 688 


27, 000 


9,938 


13,475 ! 


3,43s 


2,727 


12,425 


24,451 


23,589 


9,722 



6,592 


i 75 

4 

. 

5,600 

1,975 

9,225 

9,283 I 

83 

4 

30 

■ 

25 

35 

2,094 

88 

5 

850 

488 

1,212 

3,349 

112 

i 

4 

55 

44 

66 

638 



U OF TH! 
^ 13 P(-^ 
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Table 10. Total Coliform Concentration for Mobile Bay - 
May 1962 

Loading at Mobile River Mouth = 1/5 
TC3 j^ (D.F. = 5) 

Correction Factor for E = 500 
Kj. = 0.72 day~l 


Measured Data 


Station i 


Model 


No. of I Monthly 70 % Confidence Range Calculated 


Field 

Samnling 


Average 


3t - tSx 3c + tSj 


Result 


1,515 


18,166 
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Table 12. Total Coliform Concentration for Mobile Bay - 

July 1962 II 

Loading at Mobile River Mouth = 1/6 ij 

TC^^ (D.F. = 6) 

Correction Factor for E = 500 !? 

= 0.88 day-^ j 


Station 

No. 

Measured Data 

Model 

Calculated 

Result 

No. of 
Field 
Sampling 

Monthly 

Average 

X 

70 ^ Confidence Range 

y . - tSx 

X + tSj? 

33 . 

4 

360 

15^ 

566 

249 

36 

4 

300 

144 

456 

1,272 

59 

4 

9 

6 

12 

86 

60 

4 

35 

12 

58 

176 

61 

4 

161 

69 

101 

166 

62 


100 

50 

150 

138 

65 ! 4 

20 

9 



66 

4 

4o 

15 

65 

78 

J 

67 


33 

13 

CO 

86 

75 


120 

10 

230 

' 

29 

83 

5 

10 

9 

11 

0 

88 

5 

13 

5 

21 

20 

112 

5 

40 

20 

6C 

2 


I 
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Table 13. Total Coliform Concentration for Mobile Bay - 
August 1962 

Loading at Mobile River Mouth =1/6 
TC^^ (D.F. = 6) 

Correction Factor 'vf or E = 500 
\ = 0.90 day"^ 

















































































Figure 15. Model Calculated Averages Compared \ 
at Station No. 33. 
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Figure 17. Model Calculated Averages Compared with Actual Data of Total Coliform Concentration 
at Station No. 59. 
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Figure 19. Averages Compared with Actual Data o( Total Conform Concentration 
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Figure 21. Model Calculated Averages Co 
at Station No. 65. 
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Figure 22 . Model Calculated Averages Compared with Actua 

at Station No. 66. 







Figure 23. Model Calculated Averages Compared with 
at Station No. 67. 
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Figure 24. Model Calculated Averages Compared with Actual Data of Total Coliform Cone 
at Station No. 75. 






Figure 25. Model Calculated Averages Compared with Actual Data of Total Coliform Concentratioa 
at Station No. 83. 
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Figure 26. Model Calculated Averages Compared with Actual Data of Total Coliform Concentration 
at Station No. 88. 
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Figure 27 . Model Calculated Averages Compared with Actual Data of Total Coliform Concentration 
at Station No. 112. 
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THE INFLUENCE OF SYSTEM CHANGES ON MOBILE BAY BEHAVIOR 


In order to assess the impact that changing river flow rates, wind 
conditions, coliform bacteria, sediment loading concentration and 
water temperature have on the hydrodynamic and material transport pro- 
perties of Mobile Bay, a parametric study was conducted using the 
developed and verified mathematical models. The results of this study 
are discussed in the following sections, subdivided for clearer 
presentation of the material. 


Hydrodynamic and Salinity Models 

There are two major areas of interest for consideration in this 
section. These areas include the effects of several model parameters 
and the effect of variations in natural phenomena on the system behavior . 
Specifically, the model parameters studied are: (1) convective 

acceleration, (2) Coriolis force, (3) a non-traditional method for cal- 
culating the resistance term and (4) simulation of the salinity wedge 
effect. To weabure the influence that each of these parameters had on 
the hydrodynamic and salinity transport behavior of the bay, three levels 
of fresh water flow (12,000 cu.ft./sec., 44,000 cu.ft./sec., 245,000 
cu.ft./sec.) and two levels of wind (15 knots, 25 knots; both from the 
prevailing wind direction) were investigated. For purposes of dis- 
cussion, river flows will be held constant and comparisons made among 
the other parameters. Following comparisons at each river flow , 
observations will be made relating to the effects of variations in fresh 
water input. Tables contrasting tidal conditions, tide height curves 
predicted by the model superimposed with the actual data used in the 
verification study and salinity profiles will be presented. 


Low River Flow - The lowest river flow considered was 12,000 cubic feet 
per second for the Mobile River complex plus an additional 5000 cubic 
feet per second superimposed on Dog River. Tidal information for high 
and low water is presented in Table 15 for State Docks , Point Clear , 

Fowl River and Bon Secour River. The first column in Table 15 and all 
subsequent tables is considered to be the standard for comparison in 
that table. 

Wind direction throughout these studies was held constant at an 
angle of 200° measured clockwise from the north (southwest wind) . This 
particular direction was chosen since it approximates the prevailing 
winds in this area. Magnitudes of velocity were set at fifteen and 
twenty-five knots. It will be noted in the tables that wind from the 
southwest increases both maximum and minimum tidal amplitudes. 

Variations in salinity profiles are most significant. It can be 
noted in Figure 28 that the southwest wind pushes the fresh water toward 
the east and thus increases the salinity intrusion up the western side 
of the bay. It should be noted, also, that little happens in the 
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Table 15. Extreme Tidal Elevations (Eeet) from Mean Sea Level 
fox Low River Elow Conditions (12,000 cfs) . 


Run Number 

1 

2 

3 

4 

5 

6 

Wind Condition 

Speed, knots 

0 

15 

25 

0 

0 

0 

Direction 

SW 

SW 

SW 

SW 

SW 

SW 

Convective Acceleration 

YES 

YES 

YES 

NO 

YES 

YES 

Coriolis Force 

YES 

YES 

YES 

YES 

NO 

YES 

Resistance Calculation 
Standard 

YES 

YES 

YES 

YES 

YES 

NO 

Modified 

NO 

NO 

NO 

NO 

NO 

YES 

State Docks 

High Tide 

2.57 

2.78 

3.77 

2.60 

2.57 

2.65 

Low Tide 

-0.2d 

0.30 

1.67 

-0.42 

-0.24 

-0.50 

Point Clear 

High Tide 

2.53 

2.66 

3.21 

2.53 

2.53 

2.54 

Low Tide 

-0.27 

-0.02 

0.80 

-0.36 

• 

O 

1 

-0.41 

Fowl River 

High Tide 

2.51 

2.61 

3.05 

2.52 

2.52 

2.54 

Low Tide 

-0 . 23 

-0.01 

0.67 

-0.33 

-0.22 

-0.39 

Bon Secour River 

High Tide 

2.45 

2.51 

2.74 

2.45 

2.46 

2.50 

Low Tide 

-0.20 

-0.12 

0.29 

-0.26 

-0.17 

-0.35 
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eastern portion of Bon Secour Bay. This aspect should be kept in mind 
and observed in the other parametric studies as comparisons are made. 

The significance of this phenomenon lies in considerations of this area 
for proposed industrial developments, power generating plants or 
municipal discharges as flushing characteristics are likely to be minimal 

Referring again to Table 15 the effects of convective acceleration, 
Coriolis force and the modified resistance term can be assessed in terms 
of tidal amplitudes at high and low water. Comparing Run 1 with Run 4, 
it can be seen that considerable variations in tidal conditions exist. 

The impact on salinity profiles can be observed in a comparison of 
curves A and B in Figure 29 . 

The Coriolis force appears to have very little effect on tidal 
elevations, comparing results for Run 1 with Run 5. Surprisingly, 
however, the effect on salinity appears to be significant. This is 
especially true for the isohaline for twenty-five parts per thousand as 
can be noted by curve C in Figure 29. Since the ultimate use of the 
hydrodynamic model is for the purpose of generating data for transport 
models such as salinity, dissolved oxygen and biochemical oxygen demand, 
the Coriolis force is judged to be essential in the Mobile Bay model. 

The modified resistance term is questionable in some respects and 
promising in others. Model stability appears to be unaffected and some 
savings in computing time can be realized. However, results indicate 
that resistance has been minimized while the opposite was expected. 
Salinity effects are equally surprising. Comparing curves A and D in 
Figure 29, it will be noted that curve D for the modified resistance 
term is a closer approximation of field data than was the verification 
profile which is shown as curve A. 


Medium River Flows - An intermediate river flow for this study was 
chosen to be 44,000 cubic feet per second for the Mobile River complex 
and 5,000 cubic feet per second for Dog River. This flow regime is 
analogous to that used in the verification run for the hydrodynamic 
model . 

Wind effects on tidal amplitudes at medium river flows are very 
similar to those at low flow as was expected. Numerical comparisons can 
be made in Table 16. It should be noted that effects of the fifteen 
knot wind on tidai. amplitudes are essentially negligible for Point 
Clear, Fowl River and Bon Secour while there is significant changes at 
State Docks. Tidal effects are expected to be magnified at the latter 
location due to the convergence of land boundaries in the northern part 
of Mobile Bay. 

Variations in salinity patterns may be noted from a comparison of 
curve A with curve B and C in Figure 30. Again, it is interesting to 
note the apparent shift in fresh water flow toward the east that results 
from wind drag. The shape of the isohalines are essentially reversed 
from a generally parabolic profile pointing north for low winds to a 
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Table 16. Extreme Tidal Elevations (Feet) from Mean Sea Level 
for Medxum River Flow Conditions (44,000 cfs). 


Wind Condition 

Speed, knots 
Direction 

Convective Acceleration 

Coriolis Force 


Resistance Calculation 
Standard 
Modified 


0 

SW 


YES 

NO 


YES 

NO 


8 

9 

10 

11 

12 

15 

25 

0 

0 

0 

SW 

SW 

SW 

SW 

SW 

YES 

YES 

NO 

YES 

YES 

YES 

YES 

YES 

NO 

YES 


State Docks 


YES 

NO 


YES 
■ NO 


YES 

NO 


NO 

YES 


High Tide 
Low Tide 

Point Clear 

High Tide 
Low Tide 

Fowl River 


2.60 

2.84 

3 .87 

2.61 

2.60 

2.64 

0,05 

0.57 

1.87 

-0.12 

0.05 

-0.24 

2.50 

2.67 

3.23 

2.50 

2.51 

2.50 

-0.19 

0.01 

0.83 

-0.29 

r— 1 

o 

1 

-0.35 


High Tide 2.49 

Low Tide -0.14 


2.62 3.07 

0.02 0.72 


2.49 2.50 

-0.25 -0.15 


2.51 

-0.33 



High Tide 
Low Tide 


2.41 2.50 2.75 2.40 2.43 2.45 

-0.14 -0.11 0.31 -0.20 -0.12 -0.21 




u/ii3 
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so^^f urtaen UUomLLs for the fifteen hnot wind and app—r^tely 

centration measurements in the bay. 

Referring again are 

sLSirto'ILsrobrerved previously, ^here appears to be^a_^^^^^^ 
contribution by tbe convective acceleratxon terms a 
contribution by tbe Coriolis force term. 

Effects of tbe convective acceleration and Coriolis terms on tbe 

rin rifnrl“l!“ 

more pronounced in tbe lower portxon mriolis force 

this term is expected to have f^^^l^anges are noted in tbe 

effects are considerably more subtle. Slxgb § ^ practical 

shape of tbe curves, but tbe dxstrxbutxon of salt xs, tor a f 

purposes, tbe same. 

Results affected^byae modified 

Ilevaticns ar; lewer »ith the modified resistance term than 
those demonstrated by tbe standard metho . 

Qalirn't-v urofiles using tbe modified resistance (curve D, Figure 31) 
are nS^trsSlLr tn S-fune. 

During Aprxl of 19/3, riows reauiau t-w neolocical Survey, 

feet Der second as reported unofficially by the U. _S. Geoiogxc^ ou y 

=°"V“11rnd^UaSa“f:rthin:r‘S -It 

rStritirns’ntaf the-«aln Pass ranged f- ^rS^r^^rird^al: 
surface to thirty-two parts nor°aken “er a tidal 

limited to the area adjacent to Main P ouantitative check on the 

ri‘dS fiefah:rargfajS^^^^^^^^^^^ 

shown in Table 17 and Figures 32, 33, 34, respectxve y. 

r . tpoWi Q 17 -it can be noted that wind effects are 

slgnlflcrnrStrrlfpict^to tJdreUvatlons. Changes are somewhat more 
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Table 17 . Extreme Tidal Elevations (Feet) from Mean Sea Level 
for High River Flow Conditions (245,000 cf s) . 


Run Number 

13 

14 

15 

16 

17 

18 

Wind Condition 

Speed, knots 

0 

15 

25 

0 

0 

0 

Direction 

sw 

SW 

SW 

SW 

SW 

SW 

Convective Acceleration 

YES 

YES 

YES 

NO 

YES 

YES 

Coriolis Force 

YES 

YES 

YES 

YES 

NO 

YES 

Resistance Calculation 

Standard 

YES 

YES 

YES 

YES 

YES 


Modified 

NO 

NO 

NO 

NO 

NO 

YES 

State Docks 

High Tide 

3.60 

3.29 

4.28 

3.07 

3.05 

3.07 

Low Tide 

2.12 

2.23 

2.84 

2.04 

2.13 

1.93 

Point Clear 

High Tide 

2.55 

2.70 

3.33 

2.54 

2.57 

2.54 

Low Tide 

0.05 

0.26 

1.09 

-0.11 

0.08 

-0.17 

Fowl River 

High Tide 

2.58 

2.70 

3.21 

2.57 

2.57 

2.57 

Low Tide 

0.14 

0.31 

1.01 

0.00 

0.11 

-0.09 

Bon Secour River 

High Tide 

2.47 

2.48 

2.81 

2.43 

2.50 

2.45 

Low Tide 

-0.04 

0.07 

0.47 

-0.17 

0.00 

-0.23 
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difficult to observe in the figures for State Docks due to the depressed 
amplitudes resulting from high river flow. 

Effects of the salt wedge at low wind conditions (15 knots) may be 
observed from a comparison of curves A and B in Figure 32. Curves A and 
B represent the same data with the exception of the salinity wedge which 
was omitted in curve A. Even though significant data is not available 
for substatiation, it is believed that the profile shown as curve A is 
more representative of the actual system under high river flow conditions. 
This would indicate the need for additional data to effectively use the 
salt wedge concept in its present form at high fresh water flow con- 
ditions. Effects of variable wind may be seen from a comparison of 
curves B and C with A in Figure 33. Fresh water flows are shifted to 
the east resulting in an inversion of the salinity curves. 

Effects of the convective acceleration and Coriolis force terms on 
tidal amplitudes are analogous to those observed previously for other 
fresh water flows as can be noted in Table 17. 

Variations in salinity patterns for high river flow may be noted 
from a comparison of curve A with curve B and C in Figure 34 for con- 
vective acceleration and Coriolis force respectively. Without convective 
acceleration, isohalines near Main Pass are shifted somewhat to the east. 
Major differences are noted in isohalines for one and five parts per 
thousand. The major changes resulting from the deletion of the Coriolis 
force are observed in the isohaline for five parts per thousand as shown 
in curve C. The area that lies between the isohallnes for five and 
fifteen parts per thousand diminishes by approximately seventeen square 
miles when the Coriolis force is deleted. 

Effects of the modified resistance terms are similar to those 
previously described for low and medium river flows. Larger quantities 
of water are allowed to enter and leave the bay and this creates a tidal 
elevation greater than that observed at high tide and a tidal elevation 
less than that observed at low tide. It is expected that cross- 
sectional areas of Main Pass and Cedar Point could be modified slightly 
and possible slight variations in the depth could be made to compensate 
for the excess flows observed using the modified resistance terms. 
Variations in the salinity profile may be observed from a comparison of 
curve A with curve D in Figure 34 . The northern tip of the isohaline 
for one part per thousand is about four kilometers further north than 
the results obtained using the modified resistance method. The eastern 
portion of the five parts per thousand profile is shifted to the north 
four to six kilometers. In addition, variations are noted in the shape 
of the profiles adjacent to Main Pass. 


Summary of River Flow and Wind Influences - The effect of these system 
changes on extreme (high and low) tidal elevations at eleven locations 
(Figure 35) in the bay are shown in Table 18. In each case there is a 
pronounced influence of wind speed on tidal elevation, especially toward 
the northern bay. This is caused by the retention of water due to wind 





Table 18. Effect of Wind and River Flow on Net Velocities (FPS) 
and the Direction (Degrees) of Currents. 


River 

Flow 


(CFS) 

12,000 

1 

.077 

256.5 

2 

.018 

12.5 

3 

.038 

308. 

7_ 

4 

.052 

273.3 

5 

.077 

258.7 

6 

.035 

198.4 

7 

.152 

292.5 

8 

.207 
104 0 

Lower B 
9 

.041 
27R R 

12,000 

.051 

.039 

.072 


.055 

.119 

.034 

.123 

.205 

.047 


248.0 

305.7 

279. 

6_ 

272.1 

258.4 

202.1 

255.4 

107 3 

27 Q 7 

12,000 

, 048 

.153 

.205 


.068 

.274 

,025 

.512 

.244 

.051 


108.1 

286.3 

262,. 


263.3 

259.9 

216.9 

180.1 

117 6 

97 c; c 










1 

^ / D ,0 

44,000 

. 172 

.075 

.077 


.099 

,125 

.059 

.137 

,170 

,067 


258.6 

313.9 

288.! 

D 

y 

273.5 

264.5 

254.3 

256,9 

114 8 

252 

44,000 

. 142 

.106 

.118 


.104 

.170 

.061 

.166 

.174 

.075 


256.6 

301.1 

276,: 

l_ 

272.8 

262.2 

254.8 

227.7 

117.0 

256 1 

44,000 

.042 

.214 

,252 


.113 

.324 

.058 

.595 

.216 

.077 


216,3 

288,8 

263.2 

t- 

265.4 

258.6 

262.1 

180.3 

125.6 

257.2 


245,000 

245,000 

245,000 



256.9 

.708 

256,4 

.466 


300.0 
.493 

298.0 
.583 


272.7 

.382 

2 70.6 

.587 


273.9 

.410 

272.7 

.392 


269.7 

‘.477 

268.4 

.657 


281.7 
.326 '■ 

282.0 

.320 


252,7 I 291.5 I 264.6 


1 10 

11 

.025 

.013 

286.3 

292.6 

.076 

.043 

279.8 

308.5 


.176 

269.3 



.186 

265.7 



.556 

196.6 

.186 

214.0 

.274 

232.4 

, 120 
223.0 

.604 

,193 

.277 

.161 

193.3 

217.0 

233.7 

239.0 

.912 

.218 

.273 

.261 

182.8 

218.9 

235.7 

248.6 


300.0 


.033 

61.1 



278.4 


.007 

180.0 
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stress conditions and the decrease of bay width from 24 miles in the 
south, to 8 miles in the north. 

Similarly, significant variations in tidal elevation are observed 
at high river flows in the north with a rapid dissapation to nearly 
normal levels at the mid-bay locations. The influence of these system 
changes on current direction and speed is likewise discussed in the 
original study document (4) . 

The influence of river flow rate on bay salinity is summarized in 
Figures 28, 30, 33. Suppression of the salinity content of the intrud- 
ing Gulf waters is observed during high riverine inflow conditions. The 
bay approaches river- dominant characteristics in the upper one half of 
the bay when flows exceed ~ 150,000 cfs. Wind changes influence 
salinity distributions in the bay above 15 knot speeds. Characteristic 
shifts in fresh water flow patterns can be traced by following salinity 
profile trends from 0 to 25 knot wind conditions. Downward profiles 
gradually oscillate as the wind speed approaches 15 knots followed by a 
complete reversal in the profile at 25 knots. These shifts are directly 
related to wind stress conditions imposed by the prevailing and constant 
southwest wind investigated. 


Coliform Bacteria Distribution Model 

To assess the impact of system variable changes on coliform bacteria 
distributions in Mobile Bay, a parametric investigation was designed as 
outlined in Table 19 with corresponding results produced as listed in 
Table 20. Those parameters studied included water temperature, river 
flow rate, wind speed and direction and coliform loading at the rivers. 
Results are expressed as monthly averaged values corresponding to the 
data used for verification of the model. 


River Flow Rates - Increases in the inflow of fresh water from the river 
systems in the northern bay shift coliform bacteria group counts to the 
southeastward direction (Figure 36) . This is caused by lower retention 
times needed for the coliform bacteria to die off resulting in higher 
residual coliform concentrations in all parts of the bay. These results 
are for constant coliform loading which in most cases does not exist 
during high river flow conditions. A more realistic way of assessing^ 
the effect of changing coliform loads independent of river flow rate is 
discussed in the next section. 


Effect of Varying Coliform Loadings - Cell loading concentration of total 
coliform at the mouth of a river reflects the pathogenic pollution 
potential of the river relative to the bay. This concentration is 
attributed to waste loadings from sources such as municipal, industrial, 
and rural areas. After periods of heavy rainfall and runoff, the 
river flow rates stabilize. However, coliform loading along the river 
course usually peaks and begins decreasing at rates greater than river 




Table 19. Data used for Parametric Study runs a to r 


Hun 

kind 

River Discharge Rates 

Temperature 

Die- 
off . 
Rate 
K 

day ^ 

Speed 

knots 

e 

deg. 

Mobile 

River 

cfs 

Dog 

River 

cfs 

Tensaw 

River 

cfs 

°F 

£. 

0 

- 

24,000 

2,000 

20,000 

78.1 

0.72 

b 

i 

15 1 

1 

225 

i 

ft 

It 


ft 

ft 

! 

c i 

: 23 

223 

It 

1 

>1 

If 

1 

1 


II 

d 

15 

90 

Wt 

6> 

If 

ft 

ffl 

e 

23 

90 

ft 

If 

ft 

It 

ff 

f 

13 

315 

.. 

ft 

It 

fl 

It 

g 

23 

315 


ft 

ft 

II 

fl 

■ h 
i 

j 

0 

13 

i ' 

! 23 

223 

223 

! 

7,000 

It 

ft 

300 

tf . 1 

i 

ft 

3,000 

If 

ft 

84.2 

fl 

II 

0.90 

I* 

fl 


k 

1 

m 


n 


q. 

r 


C 

15 

25 


7.9 


225 

225 


225 


: 43 , 000 

fl 

If 

• 

3,000 

If 

tf 

i 

100,000 

ft 

II 

67.9 

ft 

If 

0.30 

i 

II 

10,000 

1,000 

9.230 

78.1 

1 0.72 

40,000 

4,000 

37,000 

ft 

II 

20,000 

2,000 

18,300 

n 

ff 

ff 

If 

If 

83.8 

0.94 

ff 

II 

fl 

68.1 

0.30 
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Table 20 . List of Figures for Parametric Study Comparisons. 


J. j 


c 

Figure 

jomparison 

among 

runs 

indicating 
the effect of 
variation of 

at 

constant 

38 

p q r 

temperature 

flow xsifos 3Xld. 

>rind 

t 

36 

nop 

river flow rate 

wind and temperature 

39 

a b c 

1 

1 

speed of wind from SW 

medium river flow 

40 

a d. e 

speed of wind from N 

medium river flow 

41 

a f g 

speed of wind from SE 

medium river flow 

1 

1 , 

42 

) 

! 

h i j 

t 



speed of wind from S<i 

I 

1 low river flow 

i 

43 

_) 

k 1 m 

speed of wind from SW 

high river flow 

i 


44 

1 

a b d f 

direction of wind at 
15 knots 

medium river flow 

■ 

45 

1 

j 

a c e g 

1 

direction of wind at 
25 knots 

medium river flow 




Figure 36. Total due to 

variatious o£ River Flow Rates. 
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flow decreases. In this discussion river flow rates, wind conditions, 
and temperature are held constant. Other parameters were set at values 
listed in Table 21. The only changes made are on the loading concentra- 
of total coliform bacteria at the mouths of Mobile River and Dog 
River. The resulting total coliform concentration profiles are shown 
in Figure 37. Comparisons are made at two concentration levels, i.e. 

70 and 1000 MPN/100 ml. Each of the shifts of the coliform concen- 
tration profile is in the order of 2 grid widths (4 km) . It is noted 
that the 70 MPN/100 ml contour shifts as many as 6 grid widths from one 
extreme to the other as 7/8 of the original total coliform bacteria is 
removed or reduced. These changes in total coliform loading are also 
more representative of conditions that might be achievable for varying 
degrees of treatment of municipal and industrial waste sources . 


Figure 38 shows the effects of changing temperatures on 
total coliform distribution. The shifts of the 100 and the 500 MPN/100 
ml total coliform concentration isolines are in the order of 2 to 4 grid 
widths (4 to 8 km.) from run to run, which can seriously affect the 
shellfish harvesting activities in the bay (10,11), especially in the 
Bon Secour areas. This simulates what can happen to the coliform dis- 
tribution in case of sharp temperature variations when all the other 
system variables, i.e, river flow rates, wind conditions, and waste 
loadings, remain unchanged. The reason for such pronounced shifts of 
concentration profiles is the change in dieoff rate constant, 
Kj., caused by temperature variation. When water temperature in the bay 
is higher, total coliform bacteria dissipate at a higher rate, and the 
coliform concentration in the bay becomes lower. When the water 
temperature is lower, is smaller, the total coliform bacteria dieoff 
at a lower rate, and the coliform concentration in the bay becomes 
higher. This effect also partly accounts for observed seasonal 
variation of total coliform concentration within Mobile Bay (12) . 


Wind Speed and Direction - Three wind speeds (0, 15, 25 knots) and three 
wind directions (north, southeast, southwest) were investigated relative 
to the effects produced in coliform bacteria distributions in the bay. 
These results are shown in Figures 39-45. In all cases wind influence 
showed the least impact on distribution patterns. This is due in part 
to the manner in which the patterns are estimated by the model (depth 
averaged) and also in part to the time frame o'>rer which calibration and 
verification events were achieved (monthly averaged data) . It is 
obvious that wind induced surface flows and mixing can influence 
bacteria distribution within the bay. However, until a stronger data 
base becomes available, one may only speculate as to the extent wind 
direction and speed alters the distribution patterns. Current model 
results are severely restricted in this area of the investigation. 
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Figure 45, Total Colifomi Concentration Profiles from Runs 
a, c, e, and g for Medium River Flow Rates; Wind 
Constantly at 25 Knots from SW, N, and SE. 
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Sediment Distribution tfodels 

Of the several models formulated for Mobile Bay, the sediment dis- 
tribution models were the most difficult because of the rather sparse 
data available for this phase of the study. As a result several different 
methods were developed to take advantage of the available data. These 
methods included correlative techniques using bathymetry for seasonal 
or longer periods of time, remote sensing photographic correlative 
techniques for within tidal and tidal cycle conditions, and more 
rigorous methods for predicting sediment ' discharge patterns from 
dredges for subsystems during tidal periods. Each of these methods 
and the corresponding results will be discussed in the following sections. 


The Effect of Seasonal Variations on Sediment Transport in Mobile Bay - 
Just as there are settling and scouring events within tidal cycles, 
there also exist seasonal variations which influence the sediment 
transportation and deposition characteristics within the bay. The 
effect of these seasonal events are studied by considering the hydro- 
dynamic behavior of the bay using mathematical modeling methods. In 
particular, for the purpose of this report., correlation of the hydro- 
dynamic and sediment transport behavior is made using a tidal cycle 
average current generated for seasonal average flow conditions. This 
technique is a convenient method of lumping variables which are diffi- 
cult to Interpret and impossible to obtain over long term periods. The 
method provides a rapid assessment of those regions more susceptible to 
high transport and/or deposition of sediment. For Mobile Bay, a value 
of 0.2 ft/sec correlates well with observed long term sediment trans- 
portation and deposition trends (i.e. a value < 0.2 ft/sec indicates a 
region of low transportation; a value > 0.2 f t/sec indicates a region 
of high transportation) . Areas which have a high transportation 
potential regardless of river flow rate are Indicated by the closed 
regions (Figure 46). These Include the Bay areas adjacent to passes and 
waters near the Mobile and Tensaw Rivers in the north 

The open areas represent regions of low sediment transportation. 
These areas include the head waters between the Mobile and Tensaw 
Rivers, the Bon Secour Bay area and regions along the western shoreline. 
The seasonal variations can be observed by following the progression of 
the high transportation potential areas from high to low river flow 
conditions. It is likely that materials deposited during low river flow 
conditions become resuspended during high river flow periods. This 
phenomenon can be traced using the hydrod 3 mamic and material transport 
model for the bay. Included in this method of analysis are allowances 
for turbulence as estimated by local dispersion coefficients. 


The Impact of Channelization on Long Term Sediment Transport in the Bay - 
In order to assess the possible long term impact that the Mobile ship 
channel has had on bay circulation and sediment transportation patterns, 
the hydrodynamic model was run under two conditions. The first set of 
conditions was derived from the 1847-1851 bay contour diagram from which 
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Figure 46 . Bathymetric Data for Mobile Bay, Alabama 











Figure 46a. Schematic Diagram Illustrating the Areas of High 
Transportation Potential at Low River Flow Conditions, 




Figure 46b. Schematic Diagram Illustrating the Areas c6 Hi oh 

irausportation Potential for Medium ll^er nrco^ditfL 




Figure 46c. Schematic Diagram Illustrating the Areas of High 
Transportation Potential for High River Flow Conditions 
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Figure 48. Sediment Deposition Trends as Reflected ly Bay Bath 3 rmetric Data for the Periods 
(A) 1852-1920 and (B) 1920-1973 (Sapp, 1975). 




Figure 49. Location Map and Transect Plan for Maintenance 
Dredging Program (Brett, 1975). 


Suspended Sediment: ConcencraCion, mg/1 



Figure 50. Comparison of Model Calculated (Dash) and Field Measured (Solid) Results for 
Maintenance Dredging Program (Brett, 1975). 
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than 0.2 ft/sec. Similar conclusions can be derived by considering the 
thickness of deposited material along the north-south sampling transect 
(Figure 50 ) . 

As the tide enters high water slack, the sediment transport was 
shifted from a north-northeasterly direction to a more easterly direction. 
Similar patterns were observed as those discussed during the flood tide 
condition. However, because of the low current velocities, the dredging 
discharge rate becomes an important source of energy for transportation 
during this period = Thus the nature of the deposition patterns was such 
that this material was deposited over a shorter distance (Figure 50 ). 

These observations are consistent with the lower velocities and shorter 
period of time that occur during the slack water condition. 

Similar patterns to those experienced during, flood tide and high 
water slack conditions were postulated for ebb tide and low water slack. 
More material will be transported over a longer distance as a consequence 
of the longer period of high current velocities in the ebb flow 
direction. During seasonal periods when river flow rates are smaller 
than the value investigated in this study, a smaller area will be 
affected as a' result of the more uniform ebb/flood tidal relationship. 
Because of the river inputs, the ebb tide condition is always greater 
than the flood conditions except during unusual periods. 


Naturally Occurring Sediment Transportation Events - The relationship of 
sediment transportation patterns to the hydrodynamic properties of 
Mobile Bay is shown by comparing model predicted velocity profiles with 
satellite photographs (Figure 51) . In the case shown the hydrodynamic 
model was run at the local conditions observed during the photographic 
mission over the bay. The resulting velocity vectors were then reduced 
by a density slicing method where the following criteria were applied: 

Category 12345 

Velocity Range 0-0.43 0.43-0.82 0.82-1.48 1*48-2.95 2.95-4.10 

ft/sec 

It should be noted that this method is highly acceptible when there 
is high sediment loads within the bay in which hydrodynamic factors con- 
stitute the primary driving force. 



/. Velocity Range 


Figure 51. Comparison of Sediment Transportation Patterns Predicted by the Hydrodynamic 
Model (A) and Skylab IV Photograph taken January 21, 1974 (B). 
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Output Format 

During the course of the previously described investigations, 
emphasis on getting the results to potential users was a primary con- 
sideration. From the outset it was planned to present model results in 
a form easily identified with the Mobile Bay system. This was achieved 
by outputing the computer results for current, salinity, coliform 
bacteria and sediment distributions in the configuration of the bay 
(see Figure 52 ) . This allows the rapid comparison of system behavior 

as a function of general location within the system. 


Report Distribution 

All of the reports were distributed to interested state, regional 
and federal agencies who could use the results provided. A listing of 
these agencies and the d-istribution for each of the interim reports is 
shown in Table 22. All of these reports were issued using contract funds 
with no expense to the agency receiving/requesting the document. On 
several occasions contacts were made by the principal investigator either 
by telephone or personal visits to the agency offices. This was 
especially the case for state and regional agencies with related water 
resource or quality interests. 


R.elated Methods of Communicating Contract Findj-ngs 

In addition to the interim report distributions, technology transfer 
was also achieved by several other methods. These included articles 
published in technical journals, presentation of results at local, 
regional or national/international meetings, conferences and workshops, 
andj seminars and interviews with various news media personnel . A listing 
of these activities are included in Table 23. In cases where results 
were presented at regional or national/international meetings, con- 
ferences and workshops, expenses were shared between contract funds and 
University funds designated to provide participation opportunities for 
faculty and research staff. 


HEIGHT FROM REFFREhCE PL»HE TO SURFJICE 
38 

37 
36 

35 • 

38 
33 
32 

31 I 

30 
29 
28 


27 





. 98 

26 





.93 

25 





.97 

28 





1.99 

23 




. 97 

.98 

22 




1 . no 

1. 01 

21 




1. 03 

1. 06 

20 



1. 07 

1. 08 

1. 09 

19 



1. 09 

1. 10 

1. 13 

19 



1. It 

1. 12 

1. 13 

17 



1. 08 

1.10 

1. IZ 

16 



1. 07 

1. 08 

1. 10 

15 

‘ 


1. 08 

1. 08 

1. 09 

18 



1. 07 

1. 07 

1. 07 

13 



1. 08 

1. 05 

1. 05 

12 

r 


1. 02 

1. 03 

I. 03 

11 

1 

. 97 

. 99 

1. 02 

1. 03 

10 

i 

1 . no 

1 . DO 

1. 01 

1.01 

9 

I. on 

. 98 

. 98 

. 99 

. 99 

8 

1 . 211 

1. Ql 

. 97 

. 96 

. 98 

7 

L 

. 97 

. 95 

. 98 

. 96 

6 




. 88 

. 90 

5 




. 76 

. 80 

8 

. 59 

. 59 

.68 

. 61 

. 68 

3 

. 57 

.611 

. 58 

. 63 

.68 

2 

.59 

. 61 

.62 

. 68 

.66 

1 

, 60 

.57 

. 59 

. 55 

. 59 


1. 37 

. 00 

.00 

.00 

.00 

. 59 


1.30 

.00 

.00 

.00 

.00 

.55 


1. 17 

.00 

.00 

. 00 

.00 

.51 


1. 05 

.00 

.00 

.00 

.00 

-SO 


.95 

.00 

.00 

.00 

.00 

.50 


. 89 

.00 

. 00 

.00 

.00 

.51 


. 89 

.00 

. 00 

.00 

.00 

. 53 


.91 

.00 

.00 

. 00 

..00 

.55 


. 93 

.00 

.00 

.00 

.00 

. 59 


. 96 

__iUU 

.00 

.00 

.00 

.65 


. 97 

.97 

ria' 

.9'8 ;5i 

-72 


. 97 

.96 

. 97 

. 97 

.96 

.95 


. 97 

. 97 

. 98 

. 98 

.98 

.98 


. 98 

.98 

. 98 

. 99 

1.00 

1.08 


.99 

.98 

. 99 

1.00 

1.06 

1.08 


1. DO 

I.Ol 

1.03 

1. 08 

1. 10 

1.11 


1.03 

1. 06 

1.12 

1. 18 

1. 18 

1. 16 

1. 19 


1. 06 

1. 18 

1.18 

1. 16 

1.16 

1.18] 



1. 13 

1.15 

1.15 

i . 1 6 

1.17 

1.20 


1. 18 ■ 

1.16 

1. 16 

1. 15 

1.16 

p.M , i 


1. 18 

1. 17 

1. 17 

1. IS 

1.15 



1. 12 

1 . 17 

1.17 

1. 16 

1.18 

1.15 


1. 11 

1.15 

1. 15 

1. 15 

1.13 

1. 13 


1. 09 

1. 12 

1. 18 

1. 13 

1.18 

1. 13 


1. 06 

1.09 

1. It 

1. 18 

1.13 

1. 13 

1. 13 

L_ 

1. US 

1. 05 

1.08 

1. it 

1.11 

1. 12 

1.11 


1. 03 

1.02 

1. 06 

1. 08 

1.09 

1.09 

1. 12 

1. 03 

1.02 

1. 05 

1. 07 

1.07 

1. 07 

1. 09 

1. 03 

1. 03 

1. 02 

1. 06 

1.06 

1.07 

1. 07 

1. 01 

1.02 

1.03 

1. 08 

1.08 

1.07 

1.03 

. 99 

1.02 

1. 06 

1.02 

1.06 

1.09 

1. 07 

. 99 

. 99 

1.08 

1. 06 

1.05 

1.07 

1.07 

. 91 

I.Ot 

.97 

I. 03 

1.05 

1-09 

1. It 

.88 .92 

. 72 r— — 

. 97 

. 98j~ 



1. 16 

.63 

. 68 

.59 

.56 

.'59 



.61 

.62 

.59 

.59 

.59 



. 55 

.57 

.57 

. 58 

.56 




Figure 52. Typical Computer Output for Mobile Bay 
Calculations (Vertical Direction 
Suppressed 2 to 1) . 
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Table 23. Transfer of Technology Developed Under NAS8-29100 via The Public Media. 


Description of Entry 

Author (s) 

Where 

Date 


Publications (7) : 

"Legal Considerations of Water Pollution in Alabama" 

April, Hill 

J. of Marine 
Science, 2 (2) 

1973 


"Adaptation of Mathematical Modeling Techniques to 
Mobile Bay for Water Quality Management" 

April, Hill 

J . of Marine 
Science, 2 (2) 

1973 


"Hydrodynamic and Material Transport Model for 
Mobile Bay, Alabama" 

April, Hill, 
& Liu 

Sjrmp . on Modeling 
Techniques 

September 1975 


"Mobile Bay Hydrodynamic and Material Transport 
Mathematical Modeling: The Importance of the 

Data Collection Base" 

"Predicting Material Transportation and 
Distribution in Mobile Bay" 

"Effects of Maintenance Dredging on Sedimentation 
in Mobile Bay, Alabama" 

April, Hill, 
Ng, Liu 

April, Brett 

Remote Sensing 
Data Users’ 
Conference 

ASCE Specialty 
Conference 

WODCON VII 
Conference 

December 1975 
January 1976 
July 1976 

.3 

CJ 
O i 

O t-3 
t?3 

"Sediment Transport and Deposition Model for 
Mobile Bay, Alabama" 

April, Brett, 
Ng 

15th Inter- 
national Con- 
ference 

July 1976 

Reports (7) : 

"Water Resources Planning for Rivers Draining 
Into Mobile Bay. Part I: Hydrodynamic and 

Salinity Models," 

April, Hill 

BER Report No . 
168-112 

January 1974 


"A Hydrod 3 mamic and Salinity Model for Mobile Bay" 

April, Hill 

BER Report No . 
169-112 

May 1974 


"Verification of the Non-conservative Species 
Model for Coliform in Mobile Bay" 

April, Liu 

BER Report No. 
174-112 

June 1974 


"Water Resources Planning for Rivers Draining 
Into Mobile Bay. Part II; Non-Conservative 
Species Transport Mobel" 

April, Liu 

BER Report No. 
185-112 

January 1975 

v£) 


Table 23 (Continued) 


Description of Entry 

Reports (Continued) : 

"Water Resources Planning for Rivers Draining 
Into Mobile Bay. Users’ Manual for the Two 
Dimensional Hydrodynamic Model" 


"Water Resources Planning for Rivers Draining 
Into Mobile Bay. Sediment Transportation 
and Distribution Patterns" 


Presentations (6) : 

"A Hydrodynamic and Salinity Model for 
Mobile Bay, Alabama" 

"A Hydrodynamic and Material Transport Model 
for Mobile Bay, Alabama" 

"Predicting Sediment Transport and Deposition 
Patterns in Mobile Bay, Alabama" 

"Predicting Material Transportation and 
Distribution in Mobile Bay" 

"Effects of Maintenance Dredging on 
Sedimentation in Mobile Bay, Alabama" 

"Sediment Transport and Deposition Model 
for Mobile Bay, Alabama" 

Seminars (8) : 

"Math Modeling Estuarine Bodies" 

"Modeling Mobile Bay - An Overall Plan" 




77th National 

April, Hill Meeting of 

AIChE June 1974 

April, Hill, Symp. on Model- 

Liu ing Techniques September 1975 

Inter . Symp . on 

April, Brett Dredging Tech. September 1975 

April, Hill, ASCE Specialty 

Ng, Liu Conference January 1976 

WODCON VII 

April, Brett Conference July 1976 

April, Brett, 15th Inter. 

Ng Conference July 1976 

April Auburn University May 1971 ^ 

April Ala. Marine 

Science Inst. July 1971 



Table 23 (Continued) 


Description of Entry 

Author (s) 

Where 

Date 

Seminars (Continued) : 

"Engineering in the Marine Environment" 

April 

Tuskegee Inst. 

April 1973 

"Mobile Bay Mathematical Modeling" 

April 

The U of A 

March 1975 


"Mobile Bay Modeling" 

April 

The U of A 

March 1975 

"Modeling an Unconventional Reaction System - 
Mobile Bay, Alabama" 

April 

U. of Tenn. 

{ 

r 

November 1975 [ 

"Material Transport in Estuarine Systems - 
The Mobile Bay Experience" 

April 

Miss. St. U. 

[: 

March 1976 li 

"Predicting Material Transport in Mobile Bay" 

April 

The U of A 

April 1976 t 

Short Courses & Conferences (2) : 

"A Model of Mobile Bay, Alabama Consortium 
of Higher Education" 

April 

ETV Series: 
Environment 

\ 

October 1973 

"Mobile Bay Hydrodynamic and Material Transport 
Mathematical Modeling: The Importance of the 

Data Collection Base" 

April 

Remote Sensing 
Data Users ' 
Conference 

1 

December 1975 i 

Interviews (2) : 
"Modeling Mobile Bay" 

April 

Capstone Comments 

August 1973 

"How Will Storms Affect Mobile's Tides?" 

April 

B'ham News 

November 1974 

Student Related Activities (4) ; 

"A Hydrodynamic and Salinity Model for Mobile Bay" 

Hill 

Ph.D. Dissert. 

1974 

"Development of Hydrodynamic and Material 
Transport Model Interface" 

Slocovich 

M.S. (non- thesis) 

1974 

"A Non-Conservative Species Transport Model for 
Mobile Bay" 

Liu 

M.S. Thesis 

VO 

1975 

i 

"Predicting Sediment Transportation in 
Mobile Bay" 

Ng 

M.S. Thesis 

1976 j 



, 




•rTnTT.rrY of the 

PAGE IS POOR 


WATER RESOURCES PLANNING 


EOR 

RIVERS DRAINING 
INTO 


MOBILE BAY 1 


CONCLUSIONS 






CONCLUSIONS 


Models for < he prediction of Lydrodynamic, salinity, coliform 
bacteria and sediment behavior have been developed for tobile Bay. These 
models allow for: (1) variability in grid size to obtain the desired 

detail in water-land boundaries, dynamic boundaries and other points of 
interest within the system, (2) the effect of the earth's rotation, 

(3) effects of rapidly changing velocities, (4) the variability in 
bottom friction, (5) resistance created by spoil banks, (6) variability 
in fresh water flow, (7) effect of wind conditions, (8) effects of the 
salinity wedge, (9) variability in the total coliform source concen- 
tration of several locations, (10) variability in water temperature as 
related to coliform bacteria dieoff rate constants, and (11) long and 
short range sediment transport and deposition trends. These paradigms 
are based on established engineering practice and available data sources, 
and constitute the necessary framework for the development of other 
models and data collection programs in support of projects related to 
water quality and water resource assessment for rivers draining into 
Mobile Bay. 

Programs for the digital computation of these models present outputs 
in a form conducive to rapid interpretation, have simple boundary con- 
ditions, and require minimum effort for utilization. Predictive 
capabilities are consistent with expectations and are adequate for 
trend analyses. These analyses may include the assessment of man's 
activities on bay behavior involving proposed dredging operations, spoil 
island construction, and pollution discharge locations to take 
advantage of optimum flushing character! tics of the bay . The models 
can also be used to study the impact that natural disturbances have on 
changes in land boundaries, wind and rainfall impacts; especially at 
levels approaching storm surge conditions . 
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ABSTRACT 


This study is the first step in the development of a comprehensive 
model for water quality measurements in Mobile Bay. Solution of the 
two-dimensional equations of change applied to Mobile Bay yields txde 
height, current patterns and salinity distribution profiles for varying 
river discharge rates (12,000-245,000 cfs) and prevailing wind conditions 
(0-25 knots, southwest). Pronounded effects are observed in the water 
movement of the Bay under the various conditions studied. 


In addition to the above results, the developed models were 
verified with field data obtained by the U.S. Army Corps of Engineers 
and others. These verification studies were necessary to evaluate the 
effect of including Coriolis force and convective acceleration in the 
model equations. A modified method for estimating the bottom friction 
term produced savings in computational time at little loss in accuracy. 
Simulation of the salt wedge effect in the Bon Secour Bay region pro- 
vided good agreement with observed isohallnes without going to the 
expense of a three dimensional model. 

The success of the hydrodynamic and salinity models obtained in this 
study has provided a good base from which BOD-DO models can be developed 
in subsequent investigations. These models will provide a baseline from 
which effects on the bay quality can be quantitatively measured. 


interim report 


on 

Contract Number NAS8-29100 


WATER RESOURCES PLANNING FOR RIVERS 
DRAINING INTO MOBILE BAY 


A HYDRODYNAMIC AND SALINITY MODEL OF MOBILE BAY 


by 

Donald 0. Hill, Research Associate 
and 

Gary C. April, Principal Investigator 


Prepared for 

National Aeronautics and ^pace A^inistration 
George C. Marshall Space 2 
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ABSTRACT 


The purposes of this research were to develop a hydrodynamic and 
salinity model for Mobile Bay, to study contributions of various terms 
in the equations developed and to evaluate the effects of certain 
natural phenomena. The equation of motion, the equation of continuity 
and the species continuity equations were reduced to two-dimensional 
form. 


Due to the complexity of the partial differential equations in- 
volved, it was necessary to effect a numerical solution utilizing finite 
differencing. Equations in this form are amenable to computer solution. 
By superimposing a grid over the system, difference equations act on 
each finite element of the grid to determine an updated estimate of 
the solution. Successive grid sweeps are made until no significant 
change in the magnitude of the parameter being studied is evidenced. 

The computerized solutions of these basic partial differential equations 
constitute the hydrodynamic and salinity models for Mobile Bay. 

Computing efficiency and perceptibility of computer output were of 
major consideration in the formulation of the Mobile Bay model. 
Specification of a constant river discharge at the northern end of the 
bay and inclusion of a marsh area above the causeway to act as a 
capacitance element simplified boundary conditions yet maintained a 
realistic result. A weir arrangement located along the two idealized 
river channels through the marsh permits drainage and overflow at 
specified tide heights. This procedure allows variation in river flow 
with only one simple change in river discharge and does not require a 
boundary change for changes in season, winds, or weather. 

Standard computational techniques for determining resistance con- 
sume considerable time, assume a significant difference in the resist- 
ance between time t and t + At and base the magnitude of the resistance 
term on the resultant of the x-y velocities . This investigation tested 
a modified computational approach for the friction terms with the 
following assumptions: the difference between t and t + At is small 

(At = 2 min.) and the resistance in the x-dlrection is dependent only 
on the x-component of velocity. While not completely satisfactory in 
predicting tidal elevations, the method shows promise as a rapid, 
accurate step when dealing with the salinity model solutions. 

The significance of convective acceleration terms and the Coriolis 
force has been considered in this research. Convective acceleration, 
in general, should be included where rapid changes in velocity are 
expected. Due to nonlinearity, an upstream differencing technique must 
be used to insure stability. The Coriolis force was shown to be 
significant as a force acting on the bay and should be retained in the 
model. 
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Development of a method showing the effects of a salt wedge in the 
salinity program made possible the predictions of salt transport without 
the excessive computing time and cost required in a three-dimensional 
model. Consideration of the salt wedge in the Bon Secour area gave an 
improved definition of isohalines at low fresh water flow conditions, but 
may introduce error at high river flows. 

Various parameters such as river flows and wind conditions were 
studied. Since the model does not make allowances for certain marsh 
areas or other land areas that become flooded at extreme tidal conditions, 
the model is not recommended for wind velocities greater than twenty- 
five knots. 

This study represents a first step in the development of a compre- 
hensive model in the Mobile Bay area. Other studies could include 
analyses of many proposed activities in the area such as: industrial 

and commercial development, municipal and urban expansion, and 
recreational development. Hopefully, implications arising from this 
study can be implemented so that the protective system for the ecology 
can keep pace with the projected growth and industrial development in 
Alabama . 
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ABSTRACT 


Total coliform group concentration data for Mobile Bay collected 
by the Alabama State Department of Health for the period from January 
1962 to August 1962 are used for the purpose of verification of the Non- 
Conservative Species Model for Mobile Bay. These coliform group con- 
centrations are obtained by analysis as described by the outlxne 
entitled "The Significance of EC Positive Organisms in Gulf Shellfish 
Growing Waters" included in the Appendix. 

The model is verified on a monthly basis, i.e. monthly average con- 
ditions are used and the model results are tabulated and compared to the 
monthly average of actual data. The standard deviations of the actual 
data are also tabulated to show the spread of the actual data. Remarks 
on whether the calculated model results fall within the ranp covered by 
the actual data are also made as a complement to the comparison between 
the calculated result and the actual data. 

The fact that the model predicts reasonable results as compared to 
the actual data is encouraging. However, the model will be further 
tested with more detailed actual data, covering a different period of 
time and will be refined in an extensive parametric study. 
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ABSTRACT 


The purpose of this research effort is to expand the mathematical 
modeling capabilities of the hydrodynamic and salinity models of Hill 
and April to include a description of non-conservative species transport 
in the Mobile Bay system. In so doing, the knowledge gained provides 
a clear insight into the effect that rivers draining into the bay have 
on water quality conditions . 

Total collform group bacteria were selected because of their 
relationship to commercial fishing ventures within bay waters . This 
item was also chosen on the basis of data availability sufficient for 
model calibration and verification. Results are presented as monthly 
average distributions corresponding to the data base used. 

In addition to the above, a parametric study was also conducted. 

In this study river flow rates, wind conditions and bay system tempera- 
tures were investigated to determine their influence on the total 
coliform concentration patterns. Of these factors temperature and 
river flow rate had a pronounced effect on the concentration profiles, 
while wind conditions showed only slight effects. Shifts in concen- 
tration profiles as much as 8 kilometers were observed in extreme cases. 

The effect of changing total coliform group loading concentrations 
at constant river flow rates and temperature was also investigated. 

As expected these loading changes had an appreciable influence on total 
coliform distribution within Mobile Bay. 

Utilization of the Non-Conservative Species Transport Model to 
predict trend behavior in the Mobile Bay system is demonstrated. Con- 
tinuing efforts to improve the data collection programs in support of 
mathematical modeling are encouraged to increase the utility and pre- 
dictive capabilities of the models. 
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effecS InfLncinrwaLrmovem^^ Sd watef quSitrin^S; 

several mathematical models based on the law^^of ^ ^ natural systems, 

«ne«u. and energy have been 'S:Lla d”irthe°'h"f 

anTi%:f „:Srj1a^L“:ri‘ SnLhed 

Environmental Applications BrancHoitLtt 

sr“:S“3£;~ 

= = FSEs= £ “ 

is basid S'^two^^dtoeMlonal’ called throughout the booklet, 

rd^-tfdret“^^^ 

dynamic model include the contlnultv x\or,ent ’"'a^ makeup the hydro- 
equations (Table 1) . ^ momentum and y-raomentum 


Tabu ^preaentatio^ ami Operational Modes of the 

Physical Models for Mobile Bay 


_ Naatt 


Eouation Form 


Continuity 

5Qx 

+ 3Qv 


5x 

ay 

Moosntun) 

apv 


x-Coceponent I 

+ go IH 


1 ■ 
1 

dx 

y-CoQponent 

spy 

+ cD 


5t 

5y 


Results 


dt 


-(R + £) 


+ Qx(2Wsitup) 


+ Qy(2Wsiu9) 


Tidal Height 


x-Componenc of 
Surface Current 


y-Cotnpoaent of 
Surface Current 


Tidal Cyclej 
Daily Avg 
[Honchly Avg 
Seasonal 

Tidal Cycle 
Daily Avg 
Monthly Avg| 
Seasonal 

Tidal Cycle 
Dally Avg 
Monthly Avg 
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transport, artificial island construction, dredging and areas subject to 
flooding are other topics which could be investigated using the mathe- 
matical modeling approach. 

To understand how the model might be. applied to these topics, this 
booklet is subdivided into four parts for the convenience of the reader. 
These are, in order: 

• Basic Concepts in Applying the Hydrodynamic Model to a Real System 

• Model Input Requirements 

• A Detailed Illustration: Application to Mobile Bay 

• Hydrodynamic Model: Program Listing 

Each section will be discussed separately, however, it is advisable that 
they be covered sequentially during the first reading to reinforce the 
basic concepts needed to understand and apply the model. 
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ABSTRACT 


The transportation and deposition of sediment play important roles 
in understanding, analysing and planning water resources and draxnage 
system programs for Mobile Bay. The material transport propertxes are 
affected by the hydrodynamic behavior of the water mass and the rivers 
draining into the bay. Channelization within the bay is also an 
important parameter that must be considered. 

The study is divided into three parts corresponding to various data 
bases identified. The first study utilizes the turbi.dity levels in the 
bay water as a natural tracer. In this case general transport patterns 
are correlated with key hydrodynamic parameters such as river flow rates 
wind conditions and bay dispersion (used as a measure of bay mixing) . 

The second part of the investigation is a localized study of sedi- 
ment and turbidity properties near the Main Pass of Mobile Bay. Data 
collected under NAS8-30810 and directed by the Marine Environmental 
Science Consortium will be analysed statistically to produce weighted 
correlations between sediment concentration and/or turbidity and key 
physical parameters. 

Thirdly, a localized investigation of sediment transportation near 
a maintenance dredge will be studied to assess the possible impact of 
man-induced disturbances. An idealized mixing reactor model wxll be 
used to calculate suspended sediment concentration as a function of 
dredge operating parameters and the bay river and tidal states. 

■ In all studies data collected from field stations or from remote 
sensing missions will be used. Discussion of the interactive ability 
of the model with a wide variety of data bases will also be presented. 
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